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Transformation  and  Self-Similarity  Properties  of 
Gamma  and  Weibull  Fragment  Size  Distributions 


Culbert  B.  Laney1 

Engility  Corp.,  8211  Terminal  Rd,  Lorton,  VA  22079  U.S.A. 


Abstract:  This  paper  describes  properties  of  two-parameter  Gamma  and  Weibull  size 
distributions,  of  the  type  commonly  used  for  liquid  and  solid  fragmentation.  Starting  with 
general  three-parameter  forms,  this  paper  systematically  derives,  categorizes,  and  compares 
logical  variations  on  two-parameter  Gamma  and  Weibull  size  distributions.  Based  on 
comparisons  with  test  data,  rare  variants  may  sometimes  perform  as  well  as,  or  better  than,  well- 
known  variants.  Gamma  distributions  have  traditionally  been  used  for  liquid  fragmentation 
while  Weibull  distributions  have  traditionally  been  used  for  solid  fragmentation.  Based  on 
comparisons  with  test  data,  Gamma  and  Weibull  distributions  appear  to  be  equally  applicable  to 
both  liquid  and  solid  fragmentation. 

Keywords:  fragment  size  distributions,  aerosol  size  distributions,  Gamma  size  distributions, 
Weibull  size  distributions,  Rosin-Rammler  size  distributions,  Nukiyama-Tanasawa  size 
distributions,  generalized  Gamma  size  distributions,  root  normal  size  distributions,  Mott-Linfoot 
size  distributions,  universal  size  distributions 

1.  Introduction 

This  paper  describes  the  transformation  and  self-similarity  properties  of  a  wide  range  of  two- 
parameter  Gaimna  and  Weibull  size  distributions,  including  those  commonly  used  for  liquid  and 
solid  fragmentation.  Transformation  properties  have  been  previously  explored  in,  e.g.,  Paloposki 
(1991).  Self-similarity  (nonnalization)  properties  have  been  previously  explored  in,  e.g.,  Lee  et. 
al.  (2004).  The  current  treatment  expands  on  these  earlier  treatments. 

Size  distributions  are  commonly  written  in  eight  different  forms.  As  defined  here,  the 
transformation  condition  requires  that  all  eight  forms  have,  at  most,  mild  (integrable) 
singularities.  It  would  not  make  physical  sense  for  a  size  distribution  to  be  well-behaved  in  one 
ordinary  fonn  but  to  experience  severe  (non-integrable)  singularities  in  another  ordinary  fonn. 

As  shown  here,  Gamma  and  Weibull  size  distributions,  as  classes,  satisfy  the  transformation 
condition  over  a  broad  range  of  parameter  space. 

Bennett  (1936)  suggested  expressing  fragment  size  distributions  in  nonnalized  fonns,  i.e.,  with 
the  independent  variable  divided  by  a  characteristic  or  average  size.  As  a  typical  outcome,  Levy 
et.  al.  (2010)  observed:  “astonishingly,  a  simple  normalization  of  the  x  axis  by  the  average 
fragment  mass  [or  volume  or  diameter]  gathers  all  the  initially  scattered  data  into  a  single  curve.” 

Friedlander  &  Wang  (1966)  noted  that  such  nonnalizations  are  simple  examples  of  self- 
preserving  or  self-similar  size  distributions.  Expanding  on  this  observation,  Spicer  &  Pratsinis 
(1996)  defined  self-similar  size  distributions  as  cases  “when  the  steady-state  size  distributions 
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scaled  by  the  average  particle  volume  [or  diameter  or  mass]  collapses  onto  a  single  size 
distribution.” 

As  defined  here,  there  are  three  self-similarity  conditions : 

1 .  When  a  size  distribution  is  expressed  in  terms  of  a  given  average  size,  it  should  actually 
obtain  that  average  size. 

2.  When  a  size  distribution  is  expressed  in  tenns  of  a  given  average  size,  and  that  average 
size  changes,  key  free  parameters  in  the  size  distribution  should  remain  the  same. 

3.  When  a  size  distribution  is  expressed  in  terms  of  a  given  average  size,  and  that  average 
size  changes,  the  size  distribution  should  stay  the  same. 

Many  previous  studies  have  concluded  that  different  materials  exposed  to  different  conditions 
may  still  experience  similar  fragmentation  outcomes.  This  suggests  the  existence  of  standard  or 
universal  size  distributions .  Well-known  examples  include  Mott  &  Linfoot  (1943),  Marshall  & 
Palmer  (1948),  and  Simmons  (1977).  In  essence,  universal  size  distributions  involve  parameters 
that  assume  a  small  number  of  discrete  values  -  as  few  as  one  or  two  -  as  opposed  to  an  infinite 
number  of  continuous  values. 

The  third  self-similarity  condition  is  an  essential  property  of  universal  size  distributions.  A 
distribution  is  hardly  “universal”  if  it  changes  when  the  arbitrary  normalization  changes.  The 
second  self-similarity  condition  is  a  desirable,  but  not  essential,  property  of  universal  size 
distributions.  Universals  are  typically  associated  with  certain  parameters.  However,  unless  a 
universal  obtains  the  second  self-similarity  condition,  these  parameters  change  when  the 
arbitrary  normalization  changes. 

2,  Size  Distributions 

Let  D  be  the  fragment  diameter  and  let  M  be  the  fragment  mass.  The  eight  most  common  ways 
of  expressing  size  distributions  are  as  follows: 

Fm  (D)  [  Fm  (M)  ]  is  the  mass  fraction  of  fragments  with  diameters  [masses]  greater  than 
or  equal  to  D  [M\. 

f M  ( D )  [  f M  ( M )  ]  is  the  mass  fraction  of  fragments  with  diameters  [masses]  in  a  range 
dD  centered  on  D  divided  by  dD  \dM  centered  on  M  divided  by  dM] 

F(D)  [  F(M)  ]  is  the  number  fraction  of  fragments  with  diameters  [masses]  greater  than 
or  equal  to  D  [M\ 

f  ( D )  [  f  (M)  ]  is  the  number  fraction  of  fragments  with  diameters  [masses]  in  a  range 
dD  centered  on  D  divided  by  dD  [dM  centered  on  M  divided  by  dM\ 

This  list  excludes  minor  variations  such  as  1  -  Fu  and  the  use  of  fragment  volume  instead  of 
fragment  mass. 
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Notice  that  Fu  is  monotone  decreasing  such  that  FM  (0)  =  1  and  FM  (co)  =  0 .  In  addition,  fM  is 
always  non-negative  such  that: 

oo 

J7m  (x)dx  =  1  (la) 

o 

Similarly,  F  is  monotone  decreasing  such  that  F( 0)  =  1  and  F(qo)  =  0 .  In  addition,  /  is  always 
non-negative  such  that: 

oo 

^f(x)dx  =  1  (lb) 

o 

In  standard  probability  theory,  F  is  called  a  complementary  cumulative  distribution  function 
(CCDF)  and  /  is  called  a  probability  density  function  (PDF). 

The  eight  fonns  defined  above  are  related  to  each  other  by  eight  equations.  The  first  of  these 
equations  is  as  follows: 


M  =  spD'" 


(2) 


where  p  is  the  density,  s  is  a  constant  shape  factor,  and  m  is  the  spatial  dimension  (1  <  m  <  3). 


For  classic  aerosols  with  nearly-spherical  droplets,  m  =  3  and  s  =  tt/ 6.  Alternatively,  Wittel  et. 
al.  (2006)  found  regular  isotropic  eggshell  fragments  have  m  =  2  and  needle-like  plate  glass 
fragments  have  m  =1.5,  which  “implies  that  fragments  have  a  self-affine  character,  meaning  that 
the  larger  they  are,  the  more  elongated  they  get.” 

The  next  two  transfonnation  equations  are  as  follows: 


Fm(M)  =  Fm(D ) 

(3) 

F(M)  =  F(D) 

(4) 

which  come  from  the  definitions  of  FM  ( D ) ,  FM  (M) ,  F(  D) ,  and  F(M) .  The  fourth 
transformation  equation  is  as  follows: 


fM  (D)  ~  Dmf(D) 


(5a) 


or  equivalently: 


fM  (M)  ~  Mf  (M) 


(5b) 


For  example,  see  Brown  (1989)  and  Brown  &  Wohletz  (1995).  The  last  four  transformation 
equations  are  well-known  and  standard: 
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‘r  dF 

Fm  ( D )  =  -  f  f  M  (x)dx ;  fM  (D)  =  -  (6) 

D  dD 

°°  Jr1 

=  f«(Jf)  =  -^| 4-  (7) 

M 

(8) 

(9) 

3.  Average  Sizes 

Consider  the  following  average  diameters: 


dF 


^(£>)  =  -J /W* ;  f' (D)  =  -  dD 


F(M)  =  -  J  f(x)dx ;  =  - 


M 


dF 

dM 


oo 


Davg=\Df(D)dD 

0 

(10) 

00 

DM^=\D.fu(D)dD 

(11) 

0 


Equation  (10)  and  (1 1)  are  known  as  the  count  mean  diameter  (CMD)  and  the  mass  mean 
diameter  (MMD),  respectively.  The  count  mean  diameter  is  also  known  as  the  arithmetic  mean 
diameter  or  the  number  mean  diameter.  Similarly: 


D’  = 


avg  oo 


AD) 

D 


dD 


D[ 


M  avg  oo 


fu(D) 

D 


dD 


(12) 


(13) 


Equation  (13)  is  known  as  the  Sauter  mean  diameter  (SMD).  Finally,  consider  the  following 
average  masses: 


Mavg=\Mf(M)dM  (14) 

0 


4 


(15) 


M„„s=\Mfu(M)dM 

0 


For  liquid  atomization,  it  is  common  to  use  the  following  ratio: 


R 


M 


D 


M 


avg 


D' 

avg 


(16) 


Notice  that  D'M  avg  emphasizes  small  fragments  while  DM  avg  emphasizes  large  fragments.  As  a 
result,  1  <  Rm  <  oo  measures  fragment  size  spread  where  RM  =  1  corresponds  to  the  least 
possible  spread  (i.e.  monodisperse)  while  RM  =  00  corresponds  to  the  greatest  possible  spread. 

As  an  empirical  observation,  different  size  distributions  are  often  almost  the  same  when  RM  is 
almost  the  same.  Thus  it  is  common  to  see  aerosol  size  distributions  specified  exclusively  in 
tenns  of  RM  ,  e.g.,  Simmons  (1977),  Wu  et.  al.  (1991),  Chou  &  Faeth  (1998),  and  Sallam  et.  al. 
(2006).  Alternatively: 


R  = 


D. 


avg 


D' 


avg 


(17) 


The  observations  made  above  about  RM  apply  equally  to  R  . 

When  transforming  size  distributions  between  different  fonns,  it  is  useful  to  have  ratios  such  as 
the  following: 


Q  = 


D 


M  avg 


D 


avg 


SM  = 


M 


M  avg 


spD 


M  avg 


s=  Mavg 


spD"' 


(18) 

(19a) 


(19b) 


Notice  that  SM  measures  skewness,  where  SM  =  1  if  the  fragment  size  distribution  is  evenly 
balanced  between  small  and  large  fragments,  e.g.,  a  unifonn  size  distribution.  Similar 
observations  apply  to  S. 
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4.  General  Three-Parameter  Form 


4.1  Introduction 

Assuming  n>  0,  Rosin  &  Rammler  (1927,  1933)  proposed  the  following  general  form  for  solid 
fragmentation: 


fM(D)~D‘  exp [-&>’]  (20) 

where  k,  n,  and  b  are  three  free  parameters.  For  a  modern  English  language  description  of  this 
distribution,  see  Stoyan  (2013).  Building  on  Rosin  &  Rammler  (1927,  1933),  Nukiyama  & 
Tanasawa  (1938a,  b)  proposed  the  following  general  fonn  for  liquid  fragmentation: 

f(D)~Dlexp[-bDB]  (21) 

where  /,  n,  and  b  are  three  free  parameters.  For  a  modem  English  language  description  of  this 
distribution,  see  Hiroyasu  (2006).  By  Equation  (5),  Equations  (20)  and  (21)  are  identical  if: 

k  -  I  +  m  (22) 

Bennett  (1936)  suggested  dividing  D  by  a  reference  diameter,  e.g.,  the  63.2%-quantile.  Using 
D  as  the  reference  diameter,  the  final  general  three-parameter  fonn  is  as  follows: 


f(D)  = 


1 

AD 


avg 


/  \ 

i 

(  \ 

n 

D 

exp 

-b 

f  D 

D  , 

D  , 

V  avs  J 

V  avs  J 

(23) 


Notice  that  Equation  (23)  appears  to  have  five  free  parameters,  namely,  A,  /,  n,  b,  and  D  . 
However,  if  Equation  (23)  obeys  Equation  (1)  and  obtains  the  correct  D  per  the  first  self¬ 
similarity  condition,  the  number  of  free  parameters  is  reduced  from  five  to  three. 

Equations  (20),  (21),  and  (23)  were  originally  derived  empirically  via  fits  to  experimental  data. 
More  recently,  Dumouchel  (2006,  2009)  showed  that  these  equations  can  be  derived  theoretically 
using  a  maximum  entropy  approach;  the  given  proof  generalizes  earlier  work,  e.g.,  Griffith 
(1943),  Li  &  Tankin  (1987),  Cousin  et.  al.  (1996). 

Besides  being  called  Rosin-Rammler  and  Nukiyama-Tanasawa  size  distributions,  Equations 
(20),  (21),  and  (23)  are  sometimes  called  generalized  Gamma  size  distributions  after  Stacy 
(1962);  see,  e.g.,  Lushnikov  (2010),  Dumouchel  (2009),  Dumouchel  et.  al.  (2012).  Less 
commonly,  Equations  (20),  (21),  and  (23)  are  called  Frechet  (1927)  size  distributions;  see,  e.g., 
Vazquez  &  Ganan-Calvo  (2010).  In  fact,  Crooks  (2010)  cataloged  over  50  different  names  for 
closely-related  distributions,  all  described  as  special  cases  of  Amoroso  (1925)  distributions.  With 
so  many  names  in  play,  Stoyan  (2013)  reasonably  suggested  that  “a  neutral  technical  name  like 
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‘powered  exponential  distribution’  might  be  more  suitable.”  However,  this  risks  adding  yet 
another  name  to  an  already  long  list. 

4.2  Transformation  Expressions 

Using  the  expressions  from  Section  2,  Equation  (23)  can  be  recast  into  eight  equivalent  fonns; 
see  Table  1.  In  Table  1,  /,  n,  b,  d,  A,  and  B  are  parameters,  in  is  the  spatial  dimension  defined  by 
Equation  (2),  T(x)  is  the  gamma  function,  and  f  (s,x)  is  the  upper  incomplete  gamma  function. 


Table  la.  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized  Gamma)  size  distributions  expressed 
in  terms  of  F  and  f. 


Table  lb.  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized  Gamma)  size  distributions  expressed 


The  expressions  in  Table  1  are  valid  if  s  >  0  and  x  >  0  in  T(y,x) .  The  first  condition  is  true  if: 

/  + 1  /  +  m  + 1 

- >0  ;  - >0 

n  n 

Equivalently: 

/  >  -1  and  n  >  0  (24a) 
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or: 


/  <  —m  - 1  and  n  <  0 


(24b) 


The  second  condition  is  true  if: 


b>0;  d  >0  (25) 

Paloposki  (1991)  first  obtained  these  results  for  the  special  case  m  —  3. 

Notice  that  Equation  (24)  changes  if  the  choice  of  the  average  fragment  size  changes.  This  is 
discussed  further  in  Section  4.3.  In  addition,  notice  that  Equation  (24)  allows  unphysical  size 
distributions  with  infinite  R  or  RM .  Section  4.4  derives  more  restrictive  conditions  that  avoid 
this. 

Notice  that  if  n  is  negative,  which  is  rare,  these  are  sometimes  called  negative  or  inverse 
generalized  Gamma  distributions,  Twomey  (1977),  Deepak  &  Box  (1982),  Kondratyev  et.  al. 
(2006).  The  best  known  negative  generalized  Gamma  (a.k.a.,  Rosin-Rammler  or  Nukiyama- 
Tanasawa)  size  distribution  is  due  to  Griffith  (1943);  see  also  Tishkoff  &  Law  (1977)  and  Grady 
&  Kipp  (1987). 

4.3  Self-Similarity  Condition  1 

The  first  self-similarity  condition  requires  that  Equation  (23)  obtain  the  correct  D  and  M  . 
Let: 


07  =  J x1+t  exp  [-  bx"  \ix  (26) 

0 


Assuming  b  >  0  and  (/  +  i  + 1)  /  n  >  0 : 


P  / + i +1' 

W,=  l  (27) 

\n\b  " 

Table  2  shows  the  parameter  settings  required  to  ensure  the  first  self-similarity  condition. 
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Table  2.  Relationships  among  parameters  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  which  ensure  the  correct  Darg  and  Mavs. 


Notice  that,  as  expected,  all  of  the  expressions  in  Table  2  are  valid  if  Equation  (24)  is  true. 
Suppose  DMavg  is  used  as  the  reference  diameter  instead  of  D  .  Then  Equation  (23)  becomes: 


f(D) 


A  D 

^M-^Mavg 


f  \ 

4/ 

(  \ 

nM 

D 

exp 

'  D 

D„ 

Du 

V  M  avg  ) 

\  M  avg  J 

(28) 


In  this  case,  Table  3  replaces  Table  2. 

Table  3.  Relationships  among  parameters  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  which  ensure  the  correct  DMavs  and  MMave. 


w  =  w 

V¥  m  vv  m+ 1 

T 

( I'M  +  m~ 1"  2 

v  nM  j 

”M  / 

A 

lM  +  m  +  \ ' 
V  nM  J 

»u 

U  CnM^m 

UM^M 

r 

Im  +2m  +  1^ 

V  nM  ) 

nM  / 

m  j 

A 

Am  +m  +  r 
l  nM  , 

nM 

m 

SM 

w2m 

Wm 

h  "m  rf  +  ^ +  ^ 

UM  1 

V  nM  J 

/r(iM+m+l] 

/  {  nM  ) 

K 

i  JmAl 1 

h  Hm  r 

UM  1 

nM 

Ou+A 

V  nM  J 

1 

nM 

lM  +m+ 1  / 

-  h  Hm  r 

|  UM  A 

\ 

l  M  +  m  +  0 

nM  J 
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Notice  that  all  the  expressions  in  Table  3  are  valid  if: 


Equivalently: 


'«±i>0  and  '«+2m  +  1>0 


lM  >  -1  and  nM  >  0 


(29a) 


lM  <  -2m  -1  and  nM  <  0 


(29b) 


In  other  words,  Equation  (29)  replaces  Equation  (24)  when  DM  replaces  D  . 

Suppose  D'  is  used  as  the  reference  diameter  instead  of  D  .  Then  Equation  (23)  becomes: 


/vro  1D  Ad 

f(D)  = -  -  exp  -b  - 

A'D'  D'  D' 

avg  \  avg  y  V  avg 


In  this  case,  Table  4  replaces  Table  2. 


Table  4.  Relationships  among  parameters  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  which  ensure  the  correct  D  'ave  and  .V/ 


V  W0  =  w, 


rf— T  M- 

{  n'J  in' 


d'  b'S' 


I'-m  +  l 


b'  ™  rf — 1  /rf  I'~m+l 


‘A  Ji'+ 1 


—  b'  n'  r 


1  jf!'  +  m  +  l 

—  b'  n  r - 


Notice  that  all  the  expressions  in  Table  4  are  valid  if: 


I'-m  +  l  .  ,  l'  +  m  +  l  . 

- - - >  0  and  - - - >  0 
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Equivalently: 


V  >  m— 1  and  ri  >  0 


(31a) 


or: 


l'  <  — m  —  1  and  n  <  0 


(31b) 


In  other  words,  Equation  (31)  replaces  Equation  (24)  when  D'  replaces  D 


avg 


Suppose  D'm  is  used  as  the  reference  diameter  instead  of  D  .  Then  Equation  (23)  becomes: 


f(D)  = 


a'  n' 

avg 


D 


n' 

avg  J 


exp 


-K 


D 


D' 

y^Mavg  J 


(32) 


In  this  case.  Table  5  replaces  Table  2. 

Table  5.  Relationships  among  parameters  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized 


Notice  that  all  the  expressions  in  Table  5  are  valid  if: 

I'm  +m  + 1 


^>0 


and 


>0 


M 


M 
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Equivalently: 


l'M  >  -1  and  nM  >  0  (33a) 

or: 

l'M  <  -m- 1  and  n'M  <  0  (33b) 

In  other  words,  Equation  (33)  replaces  Equation  (24)  when  D'Mavg  replaces  D  . 

4.4  Self-Similarity  Condition  2 

The  second  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  certain  key  parameters.  For  example,  suppose  the  invariant  parameters  are  /  and  n.  In 
other  words: 

1  =  1  =/'  =  /' 

1  1  1 M 

n  =  nM  =n'  =  n'M 

Table  6  gives  expressions  for  ratios  of  averages  for  Equation  (23),  (28),  (30),  and  (32). 

Table  6.  Ratios  of  averages  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized  Gamma)  size 
distributions  assuming  l  and  n  are  invariant. _ 
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fl  +  m  +  lV 

V  n  J 

Notice  that  the  expressions  in  Table  6  are  valid  if: 


l  A  j  l  +  m  +  2 

—  >  0  and  - >  0 
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Equivalently: 


/  >  0  and  n  >  0 


(34a) 


or: 


l  <  -m  -  2  and  n  <  0  (34b) 

Where  it  is  more  restrictive,  it  is  recommendation  that  Equation  (34)  be  used  in  place  of 
Equations  (24),  (29),  (31),  and  (33).  The  results  are  summarized  in  Table  7. 


Table  7.  Recommended  parameter  range  for  Rosin-Rammler  (a.k.a.,  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  assuming  that  /  and  n  do  not  depend  on  the  choice  of  the  average  fragment  size. 


Average  Size 

n  >  0 

n  <  0 

D aVg  ’  M avg 

/>0 

/  <  —m  -  2 

avg  ’  M m  avg 

/>0 

/  <  —2m  —  1 

D'  M' 

avg>  1V±  avg 

l>m  —  1 

/  <  —m  -  2 

D'  M' 

avg  9  iV±  M  avg 

/>0 

/  <  —m  -  2 

4.5  Self-Similarity  Condition  3 

The  third  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  the  size  distribution.  In  general,  even  when  the  size  distribution  stays  the  same,  all  the 
parameters  change;  an  iterative  technique  is  required  to  find  the  new  parameters.  However, 
suppose  that  the  parameters  /  and  n  are  invariant.  Using  Tables  2,  3,  and  6,  it  can  be  shown  that: 


0  = 


\ 

J 


1/  n 


(A 

v 


\-i/(/+i) 


M 

A 


(35) 


This,  in  turn,  implies  that  Equations  (23)  and  (28)  are  the  same.  Similarly,  using  Tables  2,  4,  and 
6,  it  can  be  shown  that: 


R  = 


f-1 


V  A  y 


l/(/+l) 


(36) 


This,  in  turn,  implies  that  Equations  (23)  and  (30)  are  the  same.  Finally,  using  Tables  2,  5,  and  6, 
it  can  be  shown  that: 


b[ 


,  -1  h, 


M 

b 


A’ 

am 


y/(/+i) 


This,  in  turn,  implies  that  Equations  (23)  and  (32)  are  the  same. 


(37) 
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5.  Alternative  Three-Parameter  Form 


5.1  Introduction 

In  Equation  (23),  replace  D  by  M  and  D  by  M  to  obtain: 


f{M) 


1 


AM 


avg 


/  X 

i 

(  \ 

n 

M 

exp 

-b 

\  M 

M 

M 

V  avs ) 

\  avg  J 

(38) 


5.2  Transformation  Condition 

Using  the  expressions  from  Section  2,  Equation  (38)  can  be  recast  into  eight  equivalent  forms; 
see  Table  8. 


The  expressions  in  Table  8  are  valid  if  s  >  0  and  x  >  0  in  r(s,x) .  The  first  condition  is  true  if: 


/+1  n  ,  1+2 
- >  0  and  - >  0 
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Equivalently: 


/  >  -1  and  n  >  0  (39a) 

or: 

/  <  -2  and  n  <  0  (39b) 

The  second  condition  is  true  if: 

b>0;d>0  (40) 

Equation  (39)  changes  if  the  choice  of  the  average  fragment  size  changes;  see  Section  5.3.  In 
addition,  Equation  (39)  allows  physical  size  distributions  with  infinite  R  or  RM;  see  Section  5.4 
for  more  restrictive  conditions  that  avoid  this. 

5.3  Self-Similarity  Condition  1 

The  first  self-similarity  condition  requires  that  Equation  (37)  obtain  the  correct  D  and  M avg  . 
Table  9  shows  the  parameter  settings  required  to  ensure  the  first  self-similarity  condition. 

Table  9.  Relationships  among  parameters  for  alternative  three-parameter  size  distributions  which  ensure  the 
correct  Dave  and  Mavs. 


Notice  that,  as  expected,  all  of  the  expressions  in  Table  9  are  valid  if  Equation  (39)  is  true. 
Suppose  Mm  is  used  as  the  reference  mass  instead  of  M  .  Then  Equation  (38)  becomes: 
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In  this  case.  Table  10  replaces  Table  9. 

Table  10.  Relationships  among  parameters  for  alternative  three-parameter  size  distributions  which  ensure 
the  correct  PMavs  and  MMave. _ 
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Notice  that  all  the  expressions  in  Table  10  are  valid  if: 


^±i>0  and  ^>0 


M 


M 


Equivalently: 


lM>- 1  and  nM  >  0  (42a) 

or: 

l M  <  -3  and  nM  <  0  (42b) 

In  other  words,  Equation  (42)  replaces  Equation  (39)  when  MMavg  replaces  M  .  Similar 
expressions  may  be  obtained  when  M'avg  or  M’m  replaces  Mavg  . 

5.4  Self-Similarity  Condition  2 

The  second  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  certain  key  parameters.  For  example,  suppose  the  invariant  parameters  are  /  and  n.  In 
other  words: 

l  =  lu 
n  =  nM 
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Then  Table  1 1  gives  expressions  for  ratios  of  averages  for  Equations  (38)  and  (41). 

Table  11.  Ratios  of  averages  for  alternative  three-parameter  size  distributions  assuming  l  and  n  are  invariant. 


0 

K  wVm¥l 

r 

f/  +  l\ 
l  n  J 

r 

^  /  + 1  /  m  +  2  ) 

v  n  J 

W,  wUm 

r 

( l  +  2\ 
l  n  J 

r 

'l  +  1/m  +  0 
l  n  ) 

R 

w  w 

yy  -l/myy \lm 

r 

/ 

r 

V 

/  —  1  /  m  + 1  ^ 

{  n  ) 

n  J 

w2 
yy  0 

T 

//  + 1 

V  n  , 

\2 

/ 

Rm 

w  w 

yy  -\lm+\yy  1/m+l 

r 

'  l  +  1/  m  +  2^ 

/ 

r 

V 

/  - 1  /  m  +  2 N 

v  n  j 

n 

r 

f!  +  2 

V  n 

y 

J 

Notice  that  the  expressions  in  Table  1 1  are  valid  if: 


/-1/m  +  l  .  ,  l  +  Hm  +  2  . 

- >  0  and  - >  0 


Equivalently: 


/  > 1/m  - 1  and  n  >  0  (43a) 

or: 

/  <-l/m-2  and  n  <  0  (43b) 

Where  it  is  more  restrictive,  it  is  recommendation  that  Equation  (43)  be  used  in  place  of 
Equations  (39)  or  (42). 

5.5  Self-Similarity  Condition  3 

The  third  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  the  size  distribution.  In  general,  even  when  the  size  distribution  stays  the  same,  all  the 
parameters  change;  an  iterative  technique  is  required  to  find  the  new  parameters.  However, 
suppose  that  the  parameters  /  and  n  are  invariant.  Then  using  Tables  9,  10,  and  1 1,  it  can  be 
shown  that  Equations  (38)  and  (41)  are  the  same.  A  similar  approach  applies  for  other  average 
sizes. 


17 


6.  Two-Parameter  Form:  Gamma  Size  Distributions 


6.1  Introduction 


Dumouchel  (2009)  and  Dumouchel  et.  al.  (2012)  note  that  the  general  form  given  in  Section  4 
has  “a  problem  of  parameter  stability  that  manifests  by  drastic  variations  of  their  values  for 
reasonable  changes  in  the  initial  conditions.”  This  and  the  following  sections  consider  possible 
two-parameter  subsets  with  more  stable  parameters.  Consider  the  following  size  distributions: 
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V  avs  J 

\  avg  J 

f(D)  = 


1 

AD 


avg 


c  \ 

b-l 

(  \ 

n 

D 

exp 

-b 

l  D 

D 

D 

V  avs  J 

\  “vg  J 

Type  III: 


fu(M)  = 


BM , 


avg 


c  \ 

b-l 

r  \ 

n 

M 

exp 

-b 

(  M 

M 

M 

V  avs  J 

\  “vg  J 

Type  IV: 


f(M)  = 


AM,. 


avg 


c  \ 

b- 1 

(  \ 

n 

M 

exp 

-b 

\  M 

M 

M 

V  avs  J 

\  avg  J 

(44) 


(45) 


(46) 


(47) 


Equations  (44)  to  (47)  will  be  referred  to  as  Gamma  size  distributions.  Type  I  and  II  Gamma  size 
distributions  are  examples  of  Rosin-Rammler  (a.k.a  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  as  discussed  in  Section  4  with: 


Type  I:  l  =  b  —  m  —  1  (i.e.,  k  =  b  —  1)  (48a) 

Type  II:  /  =  b  - 1  (48b) 

Similarly,  Type  III  and  IV  Gamma  size  distributions  are  examples  of  the  alternative  three- 
parameter  size  distributions  as  discussed  in  Section  5  with: 

Type  III:  /  =  b  —  2  (i.e.,  k  =b—\)  (49a) 

Type  IV:  l  —  b-l  (49b) 


Until  recently,  Gamma  size  distributions  were  rarely  used  for  fragmentation  or  its  reverse, 
coagulation.  In  their  first  known  appearance,  Melzak  (1953)  showed  that,  if  the  initial  size 
distribution  is  a  Type  IV  Gamma  size  distribution,  then  the  Smoluchowsk  coagulation  equation 
has  an  analytical  solution;  see  also  Friedlander  &  Wang  (1966),  Scott  (1968),  and  Lindblad 
(2005,  2007). 


18 


As  another  example,  consider  fragmentation  of  unitary  solid  or  liquid  bodies.  Based  on 
computational  results  for  Veronoi  tesselations,  Kiang  (1966)  suggested  using  Type  IV  Gamma 
size  distributions.  Since  then,  researchers  have  actively  debated  whether  two-parameter  Gamma 
size  distributions  can  be  used,  as  Kiang  conjectured,  or  whether  full  three-parameter  generalized 
Gamma  distributions  must  be  used  instead,  e.g.,  Hinde  &  Miles  (1980),  Tanemura  (1988,  2003), 
Okabe  et.  al.  (2000),  Ferenc  &  Neda  (2007). 

More  recently,  Villermaux  et.  al.  (2004)  used  Type  II  Gamma  size  distributions  to  describe 
“spray  formed  when  a  fast  gas  stream  blows  over  a  liquid  volume.”  As  a  theoretical  justification, 
they  showed  that  Gamma  size  distributions  are  analytical  solutions  of  Smoluchowski-type 
coagulation  equations,  reversing  the  role  of  input  and  output  relative  to  Melzak  (1953), 
Friedlander  &  Wang  (1966),  Scott  (1968),  and  Lindblad  (2005,  2007).  To  extend  this 
observation  from  coagulation  to  fragmentation,  they  argued  that  a  common  liquid  “fragmentation 
mechanism  ...,  somewhat  surprisingly,  consists  of  a  coalescence  process.” 

While  this  may  be  true,  there  are  at  least  two  other  possible  explanations.  First,  based  on  a 
literature  review,  Vazquez  &  Ganan-Calvo  (2010)  suggested  a  simple  “analogy  between  the 
equations  of  coalescence  and  fragmentation.”  Alternatively,  in  certain  cases,  coagulation  and 
fragmentation  may  be  inverse  processes.  Bertoin  (2006)  notes  that  “in  general,  time-reversal 
does  not  transfonn  a  fragmentation  process  into  a  coalescent  process,  and  vice-versa”  but  that  in 
selected  cases  there  may  be  “remarkable  duality  between  coagulation  and  fragmentation;”  see 
also  Dong  et.  al.  (2006).  The  properties  of  Gamma  size  distributions,  such  as  their  close 
relationship  to  Poisson-Dirichlet  distributions  as  described  in  Feller  (1971),  may  qualify  them  for 
such  duality.  If  so,  the  fact  that  Gamma  size  distributions  are  a  natural  initial  state  for 
coagulation,  as  noted  above,  means  that  they  may  also  be  a  natural  final  state  for  fragmentation. 

For  other  examples  of  Type  II  Gamma  size  distributions  see,  e.g.,  Marmottant  &  Villennaux 
(2004a,  2004b),  Bremond  &  Villermaux  (2006),  Bremond  et.  al.  (2007),  Villennaux  (2007), 
Eggers  &  Villermaux  (2008),  Villennaux  &  Bossa  (2009,  2011),  and  Lhuissier  &  Villennaux 
(2013).  In  general,  these  papers  treat  the  given  Gamma  size  distributions  as  universals.  In 
particular,  Villennaux  &  Eggers  (2008)  note  that  “universality  means  that  breakup  is  difficult  to 
control,  since  its  characteristics  are  independent  of  initial  conditions.” 

6.2  Transformation  Conditions 

Using  the  expressions  from  Section  2,  Equations  (44)  to  (47)  can  each  be  recast  into  eight 
equivalent  forms.  As  an  example,  Table  12  shows  eight  different  forms  for  Type  II  Gamma  size 
distributions.  Table  12  is  obtained  by  substituting  l  =  b  —  1  and  n  =  1  in  Table  1. 
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Table  12a.  Type  II  Gamma  distributions  for  expressed  in  terms  of  F  and /. 

r  b,  bD  /  D 

F(D)  =  .  .  “vgJ 

r  (b) 

r 

F(M)  =  - 

I '(b) 

Aim 

f\D)  = 


AD. 


avg 


r  D  A 

\Davg  j 


b- 1 


exp 


D 


D. 


avg 


f(M)  = 


mAM , 


avg 


c  \ 

b  /  m— 1 

f  \ 

1/m 

SM 

exp 

-d 

f  M 

M  , 

M  , 

V  avg  J 

\  avg  J 

Table  12b.  Type  II  Gamma  distributions  expressed  in  terms  of  FM  and  fM. 


By  Equations  (24a)  and  (48),  the  transformation  condition  for  Type  I  and  II  Gamma  size 
distributions  are  as  follows: 

Type  I:  b>  m  (50a) 

Type  II:  b  >  0  (50b) 

By  Equations  (24b),  (25),  and  (48),  b  <  0  and  n  <  0  are  specifically  excluded. 

By  Equations  (39a)  and  (49),  the  transformation  condition  for  Type  III  and  IV  Gamma  size 
distribution  are  as  follows: 


(51a) 
(51b) 

By  Equations  (38b),  (39),  and  (49),  b  <  0  and  n  <  0  are,  again,  specifically  excluded. 

6.3  Self-Similarity  Condition  1 


Type  III:  b  >  1 

Type  IV:  b>  0 


As  an  example,  let  /  =  b  —  m  —  1  in  Table  2.  Then  for  Type  I  Gamma  size  distributions,  the 
correct  D  is  obtained  if: 


b  = 


(b-m  +  p 

I/  r 

(b-m\ 

V  n  J 

/ 

l  n  J 

(52a) 


Unfortunately,  this  expression  lacks  an  analytical  solution  for  n;  solutions  must  be  found 
iteratively.  As  another  example,  let  /  =  b  —  1  in  Table  2.  Then  for  Type  II  Gamma  size 
distributions,  the  correct  D  is  obtained  if: 
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b  =  r 


(52b) 


Fortunately,  this  has  a  very  simple  analytical  solution.  In  particular,  the  relation  T(h  + 1)  =  bV(b) 
implies  n  =  1 . 


Similar  results  apply  for  Type  III  and  IV  Gamma  size  distributions;  see  Table  13  for  an  overall 
summary.  Table  13  is  derived  by  substituting  Equations  (48)  and  (49)  into  Tables  2  and  9. 


Notice  that,  as  expected,  all  of  the  expressions  in  Table  13  are  valid  if  Equations  (50)  and  (51) 
are  true. 
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Suppose  Dm  avg  is  used  as  the  reference  diameter  instead  of  D  and,  similarly,  M M  is  used 

as  the  reference  mass  instead  of  M  .  Then  Equations  (44)  to  (47)  become: 


Type  I:  fM(D)  = 


Aw  avg 


/  \ 

bM  1 

f  \ 

"m 

D 

exp 

bM 

l  D 

D., 

V  M  avg  J 

\  M  avg  J 

(53) 


Type  II:  f(D) 


A  D 

avg 


/  X 

bM  1 

(  \ 

”M 

D 

exp 

bM 

\  D 

D ,, 

D,r 

V  Mavg  J 

\  M  avg  J 

(54) 


Type  III:  fM(M) 


avg 


f  lf  \ 
M 

y^Mavg  J 


bu- 1 


exp 


-b 


f  \n» 

M 


M 


avg  J 


(55) 


Type  IV:  f(M) 


/  \6m-1 

M 


^M-^M  avg 


M avg  J 


exp 


-b 


f  v« 

M 


M 


avg  J 


(56) 


As  an  example,  let  lM  =bM  —m  —  1  in  Table  3.  Then  for  Type  I  Gamma  size  distributions,  the 
correct  DM  is  obtained  if: 


bu  =rl 


bM  +1 

-  r 

f  u  \ 

bM 

v  nM  J 

/ 

\UM  g 

(57a) 


Fortunately,  this  has  a  very  simple  analytical  solution.  In  particular,  the  relation 
r(bM+l)  =  bMT(bM)  implies  nM  =1. 


As  another  example,  let  lM  =  bM  —  1  in  Table  3.  Then  for  Type  II  Gamma  size  distributions,  the 
correct  DM  is  obtained  if: 


bM=  n 


' bM  +m  +  V 

"7r 

( i  \ 

bM  +  m 

V  nM  J 

/ 

V  nM  j 

(57b) 


Unfortunately,  this  expression  lacks  an  analytical  solution  for  nM  ;  solutions  must  be  found 
iteratively.  Similar  results  apply  for  Type  III  and  IV  Gamma  size  distributions;  see  Table  14  for 
an  overall  summary. 
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Similar  expressions  may  be  obtained  for  M'avg,  D'avg  and  M'Mavg,  D'Mavg 


6.4  Self-Similarity  Condition  2 

The  second  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  key  parameters.  Sections  4.4  and  5.4  assumed  that  /  and  n  are  invariant  while  b  is  non¬ 
invariant.  However,  for  Gamma  distributions,  /  is  a  function  of  b.  Since  b  is  non-invariant,  /  must 
also  be  non-invariant.  Thus,  in  fact,  Gamma  size  distributions  have  no  invariant  parameters. 

To  distinguish  between  the  different  parameters  sets,  A,  B,  C,  and  D  may  be  used  to  refer  to 
normalization  by  D M avg ,  Davg,  D'Mavg,  and  D'avg ,  respectively.  For  example,  Table  14  gives 

expressions  for  ratios  of  fragment  averages  for  Type  IA,  IIB,  IIIA,  and  IVB  Gamma  size 
distributions. 
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Table  15.  Ratios  of  averages  for  selected  Gamma  size  distributions. 


IA 

(  Dm  avg  ) 

IIB 

(Davg) 

1 1 1 A 

(  Dm  avg  ) 

IVB 

(Davg) 

0 

bM 

b  +  m 

bM  -1+1 lm 

b  + 1  /  m 

1 

3 

b 

bM~l 

b 

R 

bM~m 

bM~m- 1 

b 

b- 1 

(  P 

r  bM- 1+- 

V  mj 

F 

(  n 

bu~  1  — 

f  m) 

(  o 

r  b  +  - 

\  m) 

r 

l  m) 

rfc,  -i)2 

T(bf 

Rm 

bM 

bM~l 

b  +  m 

b  +  m-l 

(  n 
r  K  +- 

V  m) 

r 

(h  n 

bM 

K  m) 

(  P 

r  h  +  i  +  — 

V  mj 

r 

f  n 

b+ 1 — 

l  mj 

rfe,)2 

r(h  +  i)2 

6.5  Self-Similarity  Condition  3 

As  already  noted,  the  parameters  in  the  Gamma  size  distribution  change  when  the  normalization 
changes.  However,  the  Gamma  size  distribution  itself  may  remain  approximately  the  same.  In 
fact,  Type  IIA,  IIB,  and  IID  Gamma  size  distributions  are  nearly  identical  except  for  the  smallest 
fragments. 

Figure  1  shows  an  example  where  the  Type  IIA  curve  is  plotted  as  Q  1  f(QlD ) ,  the  Type  IIB 
curve  is  plotted  as  f(D) ,  and  the  Type  IID  curve  is  plotted  as  Rf(RD)  . 
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1E+00 


0  0.5  1  1.5  2  2.5 

D/Daug 

(b.)  Linear-linear  plane 

Figure  1.  An  example  showing  that  Gamma  size  distributions  may  be  approximately  the  same  regardless  of 
normalization.  In  particular,  the  Type  IIA  is  normalized  by  (and  ensures  the  correct  value  for)  the  mass 
mean  diameter,  the  Type  IIB  is  normalized  by  (and  ensures  the  correct  value  for)  the  count  mean  diameter, 
and  the  Type  IID  is  normalized  by  (and  ensures  the  correct  value  for)  the  average  defined  by  Equation  (12). 
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7.  Two-Parameter  Form:  Weibull  Size  Distributions 


7.1  Introduction 

Consider  the  following  size  distributions: 


Type  I: 


fuiP) 


1 


BD. 


avg 


/  \ 

n- 1 

(  \ 

n 

D 

exp 

-b 

f  D 

D  . 

D 

V  avs  J 

\  “Vg  J 

Type  II:  f(D)  = 


AD. 


avg 


r  \ 

n- 1 

r  \ 

n 

D 

exp 

-b 

f  D 

D 

D  . 

V  avs  J 

\  avg  J 

Type  III:  fu(M)  = 


BM 


avg 


c  \ 

n- 1 

(  \ 

n 

M 

exp 

-b 

\  M 

M 

M 

V  avs  J 

\  avg  J 

Type  IV:  f(M)  = 


AM 


avg 


(  .  5 

n- 1 

f  \ 

n 

M 

exp 

-b 

\  M 

M 

M 

V  avs  J 

\  “vg  J 

(58) 


(59) 


(60) 


(61) 


Equations  (58)  to  (61)  will  be  referred  to  as  Weibull  size  distributions  after  Weibull  (1939a, b). 
For  a  modern  treatment  of  Weibull  distributions,  see,  e.g.,  Rinne  (2008).  Type  I  and  II  Weibull 
size  distributions  are  examples  of  Rosin-Rammler  (a.k.a.  Nukiyama-Tanasawa  or  generalized 
Gamma)  size  distributions  as  discussed  in  Section  4  with: 

Type  I:  l  =  n  —  m— 1  (i.e.,  k  =  n  —  1)  (62a) 

Type  II:  /  =  «—  1  (62b) 

Notice  that  Equation  (62)  is  the  same  as  Equation  (48),  except  that  n  replaces  b.  Similarly,  Type 
III  and  IV  Weibull  size  distributions  are  examples  of  alternative  three-parameter  size 
distributions  as  discussed  in  Section  5  with: 

Type  III:  l  =  n  —  2  (i.e.,  k  =  n  —  1) 

Type  IV:  l  =  n  —  1 

Notice  that  Equation  (63)  is  the  same  as  Equation  (49),  except  that  n  replaces  b. 

It  can  be  shown  that  Type  I  and  III  Weibull  distributions  are  the  same,  provided  that  n  in  the 
Type  III  equals  nlm  in  the  Type  I  (and  that  the  other  parameters  are  adjusted  accordingly). 
Similarly,  Type  II  and  IV  Weibull  distributions  are  the  same,  provided  that  n  in  the  Type  IV 
equals  nlm  in  the  Type  II  (and  that  the  other  parameters  are  adjusted  accordingly). 


(63a) 

(63b) 
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Type  I  Weibull  distributions  are  commonly  known  as  Rosin-Rammler-Sperling-Bennett  (RRSB) 
distributions  after  Rosin  et.  al.  (1933)  and  Bennett  (1936).  Type  II  Weibull  distributions  are 
sometimes  known  as  generalized  Mott-Linfoot  distributions  after  Mott  &  Linfoot  (1943),  e.g., 
Mock  &  Holt  (1983),  Grady  et.  al.  (2001),  Grady  (2006),  Arnold  &  Rottenkolber  (2008). 

As  an  example,  the  meteorology  community  commonly  uses  the  Marshall-Palmer  (1948)  law  for 
raindrop  size  distributions.  The  Marshall-Palmer  law  is  a  universal  Type  II  Weibull  size 
distribution  with  n  =  1  and  m  =  3.  Villennaux  &  Bossa  (2009)  showed  that  the  Marshall-Palmer 
law  can  be  “understood  from  the  fragmentation  products  of  non-interacting,  isolated  drops.”  As 
another  example,  the  weapons  effects  community  commonly  uses  the  Mott-Linfoot  (1943) 
distributions  to  describe  metal  casing  fragments;  see  also,  e.g.,  Grady  (2006).  The  Mott-Linfoot 
distributions  are  universal  Type  II  Weibull  distributions  with  n/m  equal  to  1/2  or  1/3. 

7.2  Transformation  Condition 

Using  the  expressions  from  Section  2,  Equations  (58)  to  (61)  can  each  be  recast  into  eight 
equivalent  forms.  As  an  example,  Table  16  shows  eight  different  fonns  for  Type  II  Weibull  size 
distributions.  Table  16  is  obtained  by  substituting  /  =  n  —  1  in  Table  1.  Notice  that  the  expression 
T(l,x)  =  exp(-x)  has  been  used. 

Table  16a.  Type  II  Weibull  size  distributions  expressed  in  terms  of  F  and/. _ 


F(D)  =  exp  [~b(D/DavgY]  F(M)  =  exp  [~d  (M  /  Mavg)n,m] 


By  Equations  (24)  and  (62),  the  transformation  conditions  for  Type  I  and  II  Weibull  size 
distributions  are  as  follows: 

Type  I:  n  >  m  or  n  <  0  (64a) 

Type  II:  n  >  0  or  n<  -m  (64b) 
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By  Equations  (39)  and  (63),  the  transformation  conditions  for  Type  III  and  IV  Weibull  size 
distribution  are  as  follows: 


Type  III:  n  >  1  or  n  <  0  (65a) 

Type  IV:  n  >  0  or  n  <  -1  (65b) 

When  n  is  negative,  which  is  rare,  these  are  sometimes  called  negative  or  inverse  or  reverse 
Weibull  distributions,  e.g.,  Rinne  (2008),  Lai  (2014).  Negative  Weibull  distributions  have  mainly 
been  used  for  reliability  engineering  rather  than  fragmentation. 

7.3  Self-Similarity  Condition  1 

Table  17  shows  the  parameters  which  ensure  that  Equations  (58)  to  (61)  obtain  the  correct  D 
and  Mavg  .  Table  17  is  derived  by  substituting  Equations  (62)  and  (63)  into  Tables  2  and  9. 


Suppose  Dm  avg  is  used  as  the  reference  diameter  instead  of  D  and,  similarly,  M M  is  used 

as  the  reference  mass  instead  of  M avg  .  Assume  that  this  change  does  not  affect  the  key  parameter 
n.  Then  Equations  (58)  to  (61)  become: 
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7.4  Self-Similarity  Condition  2 

The  second  self-similarity  condition  requires  that  changing  the  average  fragment  size  does  not 
change  key  parameters.  Sections  4.4  and  5.4  assumed  that  /  and  n  were  invariant.  For  Weibull 
distributions,  /  is  a  function  of  the  invariant  n.  Thus  both  /  and  n  are,  in  fact,  invariant. 

Table  19  gives  expressions  for  ratios  of  fragment  averages.  Table  19a  is  obtained  by  substituting 
Equation  (64)  in  Table  6.  Similarly,  Table  19b  is  obtained  by  substituting  Equation  (65)  into 
Table  11. 


Table  19a.  Ratios  of  averages  for  Type  I  and  II  Weibull  distributions. 


Table  19b.  Ratios  of  averages  for  Type  III  and  IV  Weibull  distributions. 
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7.4  Self-Similarity  Condition  3 


The  proofs  given  in  Section  4.5  and  5.5  apply.  In  other  words,  when  n  is  fixed  and  the  other 
parameters  vary  as  in  Table  18,  Weibull  distributions  are  the  same  regardless  of  the  choice  of 
average  fragment  size. 

8.  Modified  Gamma  and  Weibull  Distributions 

As  defined  here,  Type  I  and  II  modified  Gamma  size  distributions  are  examples  of  Rosin- 
Rammler  (a.k.a  Nukiyama-Tanasawa  or  generalized  Gamma)  size  distributions  as  discussed  in 
Section  4  with: 


Type  I:  l  =  ib-m- 1 
Type  II:  l  =  ib—  1 


(70a) 

(70b) 


where  i  is  any  integer  greater  than  1.  Notice  that  Equation  (70)  is  the  same  as  Equation  (48), 
except  that  ib  replaces  b.  Similarly,  Type  III  and  IV  modified  Gamma  size  distributions  are 
examples  of  the  alternative  three-parameter  size  distributions  as  discussed  in  Section  5  with: 


Type  III:  l  =  ib- 2 
Type  IV:  l  =  ib- 1 


(71a) 

(71b) 


where  i  is  any  integer  greater  than  1.  Notice  that  Equation  (71)  is  the  same  as  Equation  (49), 
except  that  ib  replaces  b. 

As  defined  here,  Type  I  and  II  modified  Weibull  size  distributions  are  examples  of  Rosin- 
Rammler  (a.k.a  Nukiyama-Tanasawa  or  generalized  Gamma)  size  distributions  as  discussed  in 
Section  4  with: 


Type  I:  /  =  in  -m-\ 
Type  II:  /  =  in  —  1 


(72a) 

(72b) 


where  i  is  any  integer  greater  than  1 .  Notice  that  Equation  (72)  is  the  same  as  Equation  (62), 
except  that  in  replaces  n.  Similarly,  Type  III  and  IV  modified  Gamma  size  distributions  are 
examples  of  the  alternative  three-parameter  size  distributions  as  discussed  in  Section  5  with: 


Type  III:  /  =  in  —2 
Type  IV:  /  =  in  —1 


(73a) 

(73b) 


where  i  is  any  integer  greater  than  1.  Notice  that  Notice  that  Equation  (73)  is  the  same  as 
Equation  (63),  except  that  in  replaces  n. 
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For  example,  suppose  i  =  2.  Then  modified  Weibull  distributions  can  be  written  as  follows: 


Type  I: 


Fm(D)  = 


(  \ 

n 

(  \ 

n 

1  +  b 

D 

exp 

-b 

D 

D 

D 

\  avg  / 

\  avg  y 

Type  II: 


HD)  = 


(  \ 

n 

(  \ 

n 

1  +  b 

f  D 

exp 

-b 

l  D 

D  , 

D  , 

\  avg  J 

V  avg  J 

Type  III: 


Fm(M)  = 


(  \ 

n 

(  \ 

n 

1  +  b 

\  M 

exp 

-b 

(  M 

M 

M 

\  avg  J 

\  avg  J 

Type  IV: 


F{M)  - 


f  \ 

n 

f  \ 

n 

1  +  b 

\  M 

exp 

-b 

\  M 

M  , 

M  , 

V  avg  J 

V  avg  J 

(74) 


(75) 


(76) 


(77) 


The  expression  T(2,x)  =  (1  +  x)exp(-x)  has  been  used. 

Levy  et.  al.  (2010)  first  suggested  Type  IV  modified  Weibull  distributions.  Daphalapurkar 
(2015)  later  suggested  Type  II  modified  Weibull  distributions.  Limited  comparisons  conducted 
to  date  indicate  that  modified  Weibull  and  Gamma  distributions  are,  for  all  practical  purposes, 
identical  to  Weibull  and  Gamma  distributions.  They  are  included  here  for  the  sake  of 
completeness. 

9.  Comparisons  to  Test  Data 

This  section  compares  Gamma  and  Weibull  size  distribution  to  liquid  spray  atomization  test  data. 
Small  size  spread  data  is  taken  from  Pimentel  et.  al.  (2010),  medium  size  spread  data  is  taken 
from  Li  &  Tankin  (1987)  and  Tishkoff  (1979),  and  large  size  spread  data  is  taken  from  Tuner  & 
Moulton  (1953). 

9.1  Gamma  Size  Distributions 

Figures  2,  3,  and  4  compare  Gamma  size  distributions  to  small,  medium,  and  large  size  spread 
data,  respectively.  The  Type  II  Gamma  size  distributions  were  taken  from  Villennaux  (2007)  and 
Marmottant  &  Villennaux  (2004a,  b).  The  Type  I,  HI,  and  IV  Gamma  size  distributions  were 
obtained  by  matching  R  to  the  Type  II  Gamma  size  distributions. 

As  seen  in  these  examples,  Gamma  size  distributions  typically  match  test  data  well,  albeit  the 
agreement  sometimes  deteriorates  at  the  extremes.  With  the  parameter  choices  made  here.  Type  I 
and  II  Gamma  distributions  generally  outperform  Type  III  and  IV  Gamma  distributions. 
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Figure  2.  Gamma  size  distributions  vs.  small  size  spread  test  data  for  spray  atomization  taken  from  Pimentel 
et.  al.  (2010).  The  Type  II  parameter  b=ll  is  taken  from  Villermaux  (2007).  The  Type  I,  III,  and  IV 
parameters  are  chosen  so  that  all  size  distributions  have  7?=1.0625.  In  all  cases,  m= 3. 
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Figure  3.  Gamma  size  distributions  vs.  medium  size  spread  test  data  for  spray  atomization  taken  from  Li  & 
Tankin  (1987)  and  Tishkoff  (1979).  The  Type  II  parameter  b= 6  is  taken  from  Marmottant  &  Villermaux 
(2004a).  The  Type  I,  III,  and  IV  parameters  are  chosen  so  that  all  size  distributions  have  R=  1.2.  In  all  cases, 
m=3. 
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Figure  4.  Gamma  size  distributions  vs.  large  size  spread  test  data  for  spray  atomization  taken  from  Turner  & 
Moulton  (1953).  The  Type  II  parameter  6=2.81  is  taken  from  Marmottant  &  Villermaux  (2004b).  The  Type  I, 
III,  and  IV  parameters  are  chosen  so  that  all  size  distributions  have  7?=1.5525.  In  all  cases,  wi=3. 
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9.2  Weibull  Size  Distributions 


Figures  5,  6,  7,  and  8  compare  Weibull  size  distributions  to  small,  medium,  and  large  size  spread 
data,  respectively.  The  Type  II  Weibull  size  distributions  were  taken  from  Onose  &  Fujiwara 
(2004),  Grady  et.  al.  (2001),  and  Mott  &  Linfoot  (1943).  The  Type  I,  III,  and  IV  Weibull  size 
distributions  were  obtained  by  matching  R  to  the  Type  II  Weibull  size  distributions. 

Near  the  mean  for  small  size  spreads,  the  Weibull  distributions  skew  either  right  or  left,  while  the 
test  data  is  nearly  symmetric.  Other  than  this,  Weibull  size  distributions  offer  a  reasonably  good 
match  to  test  data.  This  is  true  despite  the  fact  that  the  chosen  Weibull  distributions  were 
originally  meant  for  solid  rather  than  liquid  fragmentation. 

Notice  that  positive  Weibull  distributions  revert  to  a  power  law,  i.e.,  they  become  linear  in  the 
log-log  plane,  for  very  small  fragments.  Similarly,  negative  Weibull  distributions  revert  to  a 
power  law,  i.e.,  they  become  linear  in  the  log-log  plane,  for  very  large  fragments. 

As  noted  earlier,  Type  I  and  III  Weibull  distributions  are  identical  to  each  other.  Similarly,  Type 
II  and  IV  Weibull  distributions  are  identical  to  each  other.  With  the  parameter  choices  made 
here,  Type  II  and  IV  Weibull  distributions  tend  to  outperform  Type  I  and  III  Weibull 
distributions. 


36 


1E+00 


IE-01 


9.  IE-02 

S' 

& 

Q 

IE -03 


IE -04 


IE -05 


- Weibull  -  Type  I  (n/m=2.33)  &  III  (n=2.33) 

- Weibull  -  Type  II  (n/m=1.5)  &  IV  (n=1.5) 

□  Pimentel  et.  al.  (2010)  -  Bosch  -  Run  1 

o  Pimentel  et.  al.  (2010)  -  Bosch  -  Run  2 

A  Pimentel  et.  al.  (2010)  -  Bosch  -  Run  3 

o  Pimentel  et.  al.  (2010)  -  Bosch  -  Run  4 

■  Pimentel  et.  al.  (2010)  -  BETE  -  Run  1 
Pimentel  et.  al.  (2010)  -  BETE  -  Run  2 
Pimentel  et.  al.  (2010)  -  BETE  -  Run  3 


Pimentel  et.  al.  (2010)  -  BETE  -  Run  4 
Pimentel  et.  al.  (2010)  -  BETE  -  Run  5 


0.01 


0.1 


1 


(a.)  Log-log  plane 


D/Davg 


—Weibull  -’ 
— Weibull  -  ‘ 
□  Pimentel  i 
Pimentel  i 
Pimentel  i 
Pimentel  i 
Pimentel  < 
Pimentel  i 
Pimentel  i 
Pimentel  i 
Pimentel  i 


?  I  (n/m=2.33)  &  III  (n=2.33) 
i  II  (n/m=1.5)  &  IV  (n=1.5) 

I.  (2010)  -  Bosch  -  Run  1 
1.(2010) -Bosch -Run  2 
I.  (2010)-  Bosch  -  Run  3 
I.  (2010) -Bosch -Run  4 
I.  (2010) -BETE -Run  1 
(2010)-  BETE- Run  2 
I.  (2010)-  BETE  -  Run  3 
I.  (2010)-  BETE  -  Run  4 
I.  (2010)-  BETE  -  Run  5 


0  0.5 

(b.)  Linear-linear  plane 


Figure  5.  Weibull  size  distributions  vs.  small  size  spread  test  data  for  spray  atomization  taken  Pimentel  et.  al. 
(2010).  All  size  distributions  have  1?=1.086  and  m=  3. 
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Figure  6.  Negative  (inverse)  Weibull  size  distributions  vs.  small  size  spread  test  data  for  spray  atomization 
taken  Pimentel  et.  al.  (2010).  The  Type  II  parameter  n/m=- 1.5  is  taken  from  Onose  &  Fujiwara  (2004).  The 
Type  I,  III,  and  IV  parameters  are  chosen  so  that  all  size  distributions  have  /f=1.086.  In  all  cases,  m=3. 
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Figure  7.  Weibull  size  distributions  vs.  medium  size  spread  test  data  for  spray  atomization  taken  from  Li  & 
Tankin  (1987)  and  Tishkoff  (1979).  The  Type  II  parameter  «/hi=0.667  is  taken  from  Grady  et.  al.  (2001).  The 
Type  I,  III,  and  IV  parameters  are  chosen  so  that  all  size  distributions  have  12=1.5708.  In  all  cases,  m=  3. 
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Figure  8.  Weibull  size  distributions  vs.  large  size  spread  test  data  for  spray  atomization  taken  from  Turner  & 
Moulton  (1953).  The  Type  II  parameter  n/m= 0.5  is  taken  from  Mott  &  Linfoot  (1943).  The  Type  I,  III,  and  IV 
parameters  are  chosen  so  that  all  size  distributions  have  7?=2.418.  In  all  cases,  m=  3. 
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9.3  Gamma  vs.  Weibull  vs.  Root  Normal  Size  Distributions 


Figures  9,  10,  and  1 1  compare  Type  II  Gamma,  Weibull,  and  root  nonnal  size  distributions  for 
small,  medium,  and  large  size  spreads,  respectively. 

In  Figure  9,  the  Gamma  size  distribution  is  taken  from  Villermaux  (2007),  the  Weibull  size 
distribution  is  taken  from  Onose  &  Fujiwara  (2004),  and  the  root  normal  size  distribution  is 
taken  from  Sallam  et.  al.  (2006);  see  also  Laney  (2015).  Notice  that  Rm=  1-05  ±  0.9%  in  all  three 
cases.  Also  notice  that  the  Gamma  and  root  normal  size  distributions  are  approximately  the 
same,  except  for  the  smallest  fragments. 

In  Figure  10,  the  Gamma  size  distribution  is  taken  from  Mannottant  &  Villermaux  (2004a),  the 
Weibull  size  distribution  is  taken  from  Grady  et.  al.  (2001),  and  the  root  normal  size  distribution 
is  taken  from  Wu  et.  al.  (1991);  see  also  Laney  (2015).  Notice  thatf%/  =  1.13  ±  0.5%  in  all  three 
cases. 

In  Figure  11,  the  Gamma  size  distribution  is  taken  from  Mannottant  &  Villermaux  (2004b),  the 
Weibull  size  distribution  is  taken  from  Mott  &  Linfoot  (1943),  and  the  root  normal  size 
distribution  is  taken  from  Empie  et.  al.  (1995,  1997);  see  also  Laney  (2015).  Notice  that 
Rm~  1-20  ±  0.7%  in  all  three  cases.  Also  notice  that  the  Weibull  and  root  normal  size 
distributions  are  approximately  the  same. 
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Figure  9.  Type  II  Gamma,  Weibull,  and  root  normal  size  distributions  vs.  small  size  spread  test  data  for  spray 
atomization  taken  from  Pimentel  et.  al.  (2010).  All  size  distributions  have  RM  ~  1.05  and  m= 3. 
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Figure  10.  Type  II  Gamma,  Weibull,  and  root  normal  size  distributions  vs.  medium  size  spread  test  data  for 
spray  atomization  taken  from  Li  &  Tankin  (1987)  and  Tishkoff  (1979).  All  size  distributions  have  RM  ~  1.13 
and  m= 3. 
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Figure  11.  Type  II  Gamma,  Weibull,  and  root  normal  size  distributions  vs.  large  size  spread  test  data  for 
spray  atomization  taken  from  Turner  &  Moulton  (1953).  All  size  distributions  have  RM  ~  1-20  and  m=3. 
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10.  Conclusions 


This  paper  has  described  a  number  of  variants  on  traditional  two-parameter  Gamma  and  Weibull 
distributions.  Each  distribution  has  been  classified  as  one  of  four  types.  Within  a  given  type,  each 
distribution  has  been  categorized  as  either  positive  or  negative,  modified  or  unmodified. 

For  Weibull  distributions,  the  traditional  choices  are  Type  I  and  Type  II,  unmodified  and 
positive;  all  other  variants  described  here  are  rare  or  new.  For  Gamma  distributions,  the 
traditional  choices  are  Type  II  and  IV,  unmodified  and  positive;  all  other  variants  described  here 
are  rare  or  new.  (Notice  that  there  are  no  negative  Gamma  distributions.)  While  initial 
comparisons  with  test  data  are  promising,  more  work  is  required  to  determine  under  what 
circumstances,  if  any,  rare  and  new  variants  equal  or  outperfonn  traditional  choices. 

Table  20  summarizes  the  transfonnation  and  self-similarity  properties  of  Gamma  and  Weibull 
distributions.  Notice  that  these  properties  are  the  same  for  all  variants.  For  additional  perspective, 
Table  20  also  includes  root  normal  size  distributions;  see  Laney  (2015). 


Table  20.  Summary  of  transformation  and  self-similarity  conditions  for  three  different  families  of  size 
distributions. 


Distribution 

Transformation 

Condition 

(No  Severe 
Singularities 
Regardless  of 
Form) 

Sell 

[-Similarity  Conditions 

1. 

(Ensures 

Correct 

Average) 

2. 

(Parameters 
Same  for 
Different 
Averages) 

3. 

(Distribution 
Same  for 
Different 
Averages) 

Gamma  (all  variants) 

V 

y 

X 

y 

Weibull  (all  variants) 

V 

y 

y 

y 

Traditional  Root  Nonnal  (Type  I) 

X 

X 

X 

y 

Improved  Root  Normal  (Type  II) 

V 

y 

X 

y 

Gamma  and  root  normal  size  distributions  are  traditionally  used  for  liquid  fragmentation  while 
Weibull  distributions  are  traditionally  used  for  solid  fragmentation.  However,  based  on  the 
comparisons  given  here,  they  both  appear  to  be  equally  applicable  to  liquid  and  solid 
fragmentation.  This  is  not  surprising  given  that  three-parameter  Rosin-Rammler  size 
distributions  -  which  underlie  Gamma  and  Weibull  size  distributions  -  are  commonly  used  for 
both  liquid  and  solid  fragmentation. 

Except  for  power  laws,  most  if  not  all  universal  size  distributions  suggested  in  the  research 
literature  belong  to  one  of  the  three  families  listed  in  Table  20.  Of  the  families  listed  in  Table  20, 
only  Weibull  size  distributions  obtain  the  second  self-similarity  condition.  This  means  that  the 
exact  parameter  choice  is  potentially  meaningful  for  those  universals,  such  as  Mott-T .infoot  and 
Marshall-Palmer,  expressed  as  Weibull  size  distributions  -  but  not  for  those  expressed  as 
Gamma  or  root  normal  size  distributions. 
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